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CN We investigate the conditional vorticity budget of fully developed three-dimensional homogeneous 

^ isotropic turbulence with respect to coherent and incoherent flow contributions. The Coherent Vor- 

, ^ ticity Extraction based on orthogonal wavelets allows to decompose the vorticity field into coherent 

_^ and incoherent contributions, of which the latter are noise-like. The impact of the vortex structures 

observed in fully developed turbulence on statistical balance equations is quantified considering 
the conditional vorticity budget. The connection between the basic structures present in the flow 
and their statistical implications is thereby assessed. The results are compared to those obtained 
I— I for large- and small-scale contributions using a Fourier decomposition, which reveals pronounced 

G differences. 
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I. INTRODUCTION 



The problem of turbulence remains a paradigm for non-equilibrium statistical mechanics. The challenge 

comes from the spatio-temporal complexity which calls for a statistical description of the phenomenon, 

ri ' moreover the presence of coherent structures induces strong statistical correlations. For example, already 

^ . the single-point vorticity statistics displays a highly non-Gaussian shape, which indicates pronounced spatial 

~~' correlations of the vorticity field. Visualizations of vorticity in fully developed turbulent flows show that 

these correlations become manifest in the form of slender vortex tubes, which form a complex entangled 

global structur#M^. Hence, one of the most interesting problems in turbulence research is to understand the 

w^ relation between the coherent structures and their implications for statistical properties of the flow. 

^\ In this context it is particularly interesting to study dynamical rather than purely kinematic statistical 

^^ relations. Maybe the most fundamental dynamical balance equation related to vorticity is the balance of 

\l enstrophy production and dissipation. Deriving an equation for the vorticity probability density function 

(T^ (PDF) is even more informative and allows to study the conditional budget of enstrophy production and 

^^ dissipation, or equivalently the balance of conditional vortex stretching and diffusion, where the ordinary 

04 budget equation is contained as a special case. This budget equation was introduced by Novikov and has 

\ I been studied in a number of publication*^"^!. The conditional vorticity budget allows to quantify vortex 

^ stretching and vorticity diffusion as a function of vorticity magnitude and hence to statistically discriminate 

J^ strong vorticity regions in the flow from weak ones. One, however, would like to go one step beyond and 

rN disentangle the influence of coherent and incoherent flow contributions on these s tatist ical quantities, which 

^ is possible with the help of the orthogonal wavelet decomposition. Farge et al.l^Iir proposed a method, 

called Coherent Vorticity Extraction (CVE), to extract the coherent structures out of turbulent flows. This 

technique is based on a denoising of vorticity in wavelet space. It was shown that the CVE method is more 

efficient than Fourier filtering^^ and that fewer wavelet coefficients are necessary to reconstruct the coherent 

structures with increasing Reynolds numbepl^lj which means that the CVE method becomes more attractive 

as the flow becomes more intermittent. This method was applied to study the vortical structures in sheared 

and rotating turbulencel^ and mixing layerglS. For all investigated flows it was shown that the coherent 

vortices are well represented with few wavelet coefficients and the statistics of the remaining background flow 

exhibits more Gaussian- like behavior. Moreover, the Coherent Vorticity Simulation (CVS), which models 
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turbulent flows by considering only the time evolution of the coherent contribution w hile n eglecting the 
incoherent one to model turbulent dissipation, was introduced by Farge and co-worker j^^^^^ l The wavelet 
and Fourier nonlinear filtering methods were compared recently^- . 

The aim of the present article is twofold. First, we combine the statistical analysis of the vorticity field 
with the Coherent Vorticity Extraction in order to obtain new insights on the coherent and incoherent 
contributions to the statistics. This gives a characterization of the statistical impact of coherent structures. 
Second, the results also serve as a benchmark to characterize the performance of the wavelet decomposition. 
For example, aiming at coherent vorticity simulations, it is especially desirable that the coherent contri- 
butions are driving the nonlinear dynamics. One of the most simple checks consequently is to study their 
contributions to dynamical budget equations such as the enstrophy budget. For comparison we include an 
analysis of low and high pass Fourier-filtered vorticity fields using the same number of degrees of freedom. 
This is an interesting investigation in its own right as it yields a characterization of large- and small-scale 
contributions and their interaction. 

The remainder of this article is structured as follows. We first review the theoretical background for 
the conditional vorticity budget, the orthogonal wavelet decomposition and the principle of the Coherent 
Vorticity Extraction. Then, after summarizing some technical details on the direct numerical simulations 
performed for this work, we will present and discuss the numerically obtained results, before we conclude. 

II. PDF EQUATION AND CONDITIONAL VORTICITY BUDGET 

The dynamics of incompressible flows can be described in terms of the vorticity field ti;(a;, t) — W x u{x, t), 
defined as the curl of the velocity field. Its evolution equation takes the form 

d 

—uj + u ■ \Iuj = So; + i/Ao; + V x F, (1) 

at 

where S(a;,f) — ^ \yu[x,t) + (VM(a;,f))^] denotes the rate-of-strain tensor, v denotes the kinematic vis- 
cosity, and F[x, t) represents an external large-scale forcing applied to the flow in order to maintain a 
statistically stationary flow. As we are dealing with incompressible flows, we additionally have 

V • M = 0. (2) 

If one is now interested in the single-point statistics of the vorticity, a comprehensive characterization can be 
obtained by studying the evolution equation of the vorticity probability density function. The vorticity PDF 
can be introduced as an ensemble average ((•)) over the fine-grained PDF (the delta distribution) according 
tcP 

f[n-x,t)^{5{uj{x,t)-n)), (3) 

by which we have introduced the sample space variable Jl. In general this PDF will be a function of space 
and time, but for statistically homogeneous and stationary turbulence the PDF becomes independent of 
these variables. For isotropic turbulence the PDF only depends on the magnitude Vl of the vorticity, such 
that the PDF of the vorticity vector, /(il), can be expressed in terms of the PDF of the magnitude of the 
vorticity, /(^i), according to 

!{n) = 47rr!2/(fl). (4) 

Consequently the single-point statistics of statistically stationary homogeneous isotropic turbulence is fully 
characterized by the PDF of the vorticity magnitude only. 

Regarding the evolution of this quantity, standard PDF methods allow to derive an evolution equation 
for the single-point vorticity PDF of homogeneous turbulence which takes the forrrP^JtSil 

^/(f2; t) = -Vn ■ [(So; + //Au; + V x F\n)f{n; t)] . (5) 

The temporal change of the vorticity PDF is accordingly given by the divergence of the conditionally 
averaged right-hand side of equation ([T]) times the vorticity PDF, such that this equation takes the form of 
a continuity equation (or conservation law) for the probability density. That means, the closure problem of 
turbulence in this formulation appears in terms of the unknown conditional averages of the vortex stretching 



term, the diffusive term and tfie externai forcing. Once these functions are known, equation (Is]) is closed. 
For isotropic turbulence these terms can be simplified further as they have to take the form 

(Sa;|rj) = s(r2,i)n s{n,t) = {Q-SLj\n), (6) 

{iyAu\fi) = d{n,t)n d{n,t) = {i^Qi-Au}\n), (7) 

(V X F\n) =e{n,t)^ e(f7,t) = ((D- (V X F)|r2), (8) 

in order to maintain invariance under arbitrary rotations (and reflections). This shows that these terms 
can be characterized by scalar- valued functions depending only on the magnitude of vorticity, which can be 
obtained by projecting the vectorial conditional averages on the direction of the vorticity, a) ~ (j^/oj, where 
uj = \\<j^\\- For a more detailed account on the exploitation of statistical symmetries we would like to refer the 
reader to Wilczek et al.^^. If we insert these expressions together with relation Q into the PDF equation 
(Is]) , we obtain the transport equation for the PDF of the magnitude of vorticity 

g-Jin;t) = -—[s{n,t) + d{n,t) + ein,t)]fin;t), (9) 

by which the problem becomes eventually one-dimensional; the temporal evolution of the vorticity PDF may 
fully be characterized by knowing the functions s, d and e, which are related to the vortex stretching term, 
the diffusive term and the external forcing term, respectively. When we now consider statistically stationary 
turbulence, this equation implies 

s{Q) + d{n) + e{n) = (10) 

as the probability current for this type of statio nary o ne-dimensional problem has to vanish. Furthermore, 
it has been shown in a number of publicationa^EEl] that the conditional vortex stretching and diffusive 
term balance at sufficiently high Reynolds numbers. Compared to those terms, the external forcing has a 
negligible effect, such that we obtain the approximation 

s{n) + d{n)^o, (11) 

or equivalently 

[s{n) + d{n)] fl = (Sa;|J7) + {iyAu\fl) w 0. (12) 

This central relation states that vortex stretching and vorticity diffusion tend to cancel for a fixed magnitude 
of vorticity on the statistical average. This balance has been extensively discussed by Novikov and co- 
workera^iSl. Recently its relation to the shape and evolution of the vorticity PDF was investigatecpl]. 

The conditional balance is much more informative than the ordinary enstrophy balance as we, for example, 
can discuss the results as a function of vorticity magnitude highlighting possible correlations. Of course, the 
average enstrophy balance (discussed, e.g., by Tennekes and Lumlcy^"^) is obtained from equation (12) by 
multiplying with ftf{fl) and integrating out 17, 

(w • Suj) + {lyu ■ Aw) ^ ft- [(Sa;|n) + {iyAuj\n)] /{ft) dn (13) 

/"OO 

= / n[s{n) + d{n)]f{n)dnf=iO. (i4) 

Jo 

Higher-order moment relations can also be obtained in the same manner. The fact that conditional averaging 
allows to study the vorticity budget as a function of the vorticity eventually permits to discriminate the 
statistics of strong vorticity events from weak vorticity events as demonstrated later on. 

If we now want to establish a connection between this conditional balance and the coherent vorticity 
structures present in the flow, we have to discriminate coherent from incoherent contributions to the vorticity 
fleld. The challenge in this context is to defl ne wh at exactly a coherent structure is, and many different 
approaches have been introduced in recent years^SESl. A conceptually different way to determine the coherent 
contributions is to analyze the vorticity fleld in terms of CVE introduced by Farge and co-workers^SEiESl 
which is based on a denoising approach. 

To this end, we decompose the vorticity fleld into coherent and incoherent contributions according to 

uj{x,t) =uj''{x,t)+uj'{x,t) (15) 



using an orthogonal wavelet decomposition as described in section III By this decomposition we obtain for 
the conditional diffusive term 

{vAuj\n) = (lyAuj^ln) + {iyAu;'\n) = [d%n) + d\n)] n, (16) 

such that we get two separate contributions from the coherent and incoherent parts of all fields of the 
ensemble. This decomposition, for instance, allows to quantify how much enstrophy dissipation is contained 
in the two terms. 

The nonlinear vortex stretching term turns out to be more complicated as the rate-of-strain tensor con- 
tains both coherent and incoherent contributions. This can be seen by noting that we can calculate the 
coherent and incoherent velocity, u'^{x,t) and u^{x,t), from the decomposed vorticity field via Biot-Savart's 
law and subsequently the coherent and incoherent rate-of-strain tensor S'^{x,t) and S^{x,t). Hence, the 
vortex stretching term may be split up into four terms containing all possible combinations of coherent and 
incoherent parts of the vorticity and the rate-of-strain tensor, 

(Sw|n) = (S^uj^ln) + (S^ij'ln) + (S'w'lfl) + (S''a;^|n) (17) 

= [s'"'{n) + s"(f]) + s'%n) + s"(f))] f2. (18) 

The interesting fact now is that the different functions characterize the interaction of coherent and incoherent 
contributions of the rate-of-strain field with coherent and incoherent contributions of the vorticity field. This 
will, for instance, eventually allow to quantify if vortex stretching is mainly caused by coherent structures 
or not. 

It is also useful to compare the results obtained by wavelet filtering to another more classical filtering 
technique. To this end we make use of Fourier decomposition into ideal low pass filtered and high pass 
filtered contributions of the vorticity field according to 

i^{x,t) =iv\x,t)+u}'\x,t). (19) 

The cutoff wavenumber is determined such that the number of Fourier modes representing the low pass 
filtered vorticity field approximately corresponds to the number of wavelet coefficients representing the 



coherent vorticity field. More details on the wavelet filtering technique will be given in section III 



III. COHERENT VORTICITY EXTRACTION 

In the following, the notation for the orthogonal wavelet decomposition of a three-dimensional vector- 
valued field is presented, and the main ideas of CVE are explained. For more details, the reader is referred 
to the original paper by Farge et al.^^. 

The vector field vorticity is considered, u;{x) = (a;i,a;2,a;3), x — (xi, 2:2, 2:3) € [0, 27r]'^ given at 
resolution N = 2'^"', where J is the number of octaves in each space direction. The decomposition of uj into 
an orthogonal wavelet series yields 

uj{x) = Y,^x^xix), (20) 

XeA 

where the multi-index A = {j,ix,'iy,'i-z, fJ-) denotes the scale index j, the position index i ~ (ia;i *j/: *z): and 
the seven directions /i = 1, 2, ..., 7 of the wavelets, respectively. The corresponding index set A is given by 

A = {A = {j,ix,iy,iz,n); j = 0, ..., J- 1; ix,iy,iz = 0, ...,2'' - 1; and n = 1,...,7}. (21) 

The orthogonality of the chosen wavelets implies that the coefficients are given by a; ^ = J,q 27ri3 '^{^)''l'x{x) dx. 
The coefficients measure fluctuations of a; around scale 2~^ and around position 27ri/2-' in one of the seven 
possible directions fi. The fast wavelet transform, which has linear complexity, is used to compute efficiently 
the A^ — 1 wavelet coefficients u)x (and the mean value which vanishes in the present case) from the N grid 
point values of u). The Coiflet 30 wavele lP^ is chosen in this study; it has 10 vanishing moments and is well 
adapted to our case. In past publication d^^ l ^^ l also Coiflet 12 wavelets have been used. It has, however, been 
checked that statistical results are robust with respect to this particular choice. The main steps of the CVE 
method and our analysis are summarized in the following: 



N Rx Urms u L T 



Table I. Simulation parameters. Number of grid points A'^, Reynolds number based on the Taylor micro-scale Rx, root- 
mean-square velocity iirms, kinematic viscosity v, integral length scale L, large-eddy turnover time T, Kolmogorov 
length scale rj, and Kolmogorov time scale r^. 

i) The wavelet coefficients of the vorticity field, ujx, are obtained by computing the fast wavelet transform 
of each component of the vorticity vector. 

ii) Then we threshold the wavelet coefficients ujx, and the field can thus be split into two contributions. 
The wavelet coefficients of the coherent part are defined as 



if \\ux\\>s^^2aj\nN ^^^^ 

else 

and those of the incoherent part as the remainder. Here u/ = y/ {{uj'^)'^) f / 3 is the standard deviation 
of the individual incoherent vorticity field ((•)/ denotes the spatial average over a single field) and N 
the total number of grid points. The threshold e is motivated as an optimal method for denoising 
signals of inhomogeneous regularity^Zl and is sometimes called 'universal' as it does not depend on the 
signal, but only on the sampling size and on the variance of the noise. However, as the variance of 
the noise, which equals two-third the enstrophy of the incoherent field, is unknown, a first threshold is 
calculated from the total field and the thresholding is applied. In the next step, the total field is then 
split with a threshold calculated from the incoherent field. An iterative algorithm was developed to 
obtain the optimal threshold by Azzalini et al.l^. In the present work we have chosen to perform one 
iteration in order to privilege a good compression rate rather than a perfectly denoised contribution. 

iii) The fast inverse wavelet transform is applied to reconstruct the coherent vorticity a;^. The incoherent 
vorticity u)^ is obtained by subtracting a;^ pointwise, such as w' = a; — u)*^. Since by construction the 
two fields w'^ and u;' are orthogonal, the separation of the total enstrophy Z — (a;^)/2 into Z = Z'^ + Z^ 
is ensured because the interaction term J,„ 2^13 u)'^{x) ■ uj^{x) dx vanishes. 

iv) Biot-Savart's law u — —V x (V^'^u;) is used to obtain the corresponding total, coherent and incoherent 
velocities, where again we have it = tt"^ + it*. However, as the Biot-Savart operator is not diagonal in 
wavelet space, the decomposition of the turbulent kinetic energy E = (m^)/2 \s E = E"^' + E"^ + €e 
where eg ^ but smalP^. 



v) From the decomposed velocity field we then can calculate the total rate-of-strain tensor S 
\ [Vm + (Vm)-'"] , as well as the coherent and incoherent contributions according to S"^ = ^ [Vit"^ + (Vit^) 



LL j I , CIB well CtO UilC UUllCiCilU CliiU liH.Uiit;i CliL (^UiiLiiUUUiUiia at_.»^UH.liii^ LU kj — rj I V u, T \\/ iJi' } 



and S' = i [Vm' + [Vu'Y] , respectively. 



vi) Finally, also the Laplacian of total, as well as the coherent and incoherent contributions of the vorticity 
field is calculated. 



IV. DIRECT NUMERICAL SIMULATION 

The presented flow is generated with a standard, dealiased Fourier pseudospectral code^SESl for the vortic- 
ity equation. The integration domain is a triply periodic box of side- length 27r. The time stepping scheme 
is a third-order Runge-Kutta scheme-^. For the statistically stationary simulations a large-scale forcing is 
applied to the flow. Here, care has to be taken in order to fulfill the statistical symmetries we make use 
of in our theoretical framework^^. After numerous tests, we chose a large-scale forcing which conserves the 
kinetic energy of the flow by amplifying the magnitude of Fourier modes in a wavenumber band and letting 
their phases evolve freely. This forcing has been found to deliver satisfactory results concerning the statis- 
tical symmetries. The numerical results presented in the following are obtained from a simulation whose 
simulation parameters are listed in table ll] 
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Figure 1. Volume visualization of a wavelet-decomposed (upper panels) and Fourier-decomposed (lower panels) vor- 
ticity field (shown are the magnitudes). The upper left shows the total field, the pictures in the middle show the 
coherent/low pass filtered contribution, whereas the right pictures show the incoherent/high pass filtered contribu- 
tions. The lower left picture indicates the different color scales used (the standard deviation of the ensemble of 
total fields is cr = y/{u}^)/3 = 5.87). The visualizations have been produced with the free software package VAPOR 
(www . vapor . ucar . edu) . 

For the numerical evaluation we average over twenty realizations of the vorticity field to ensure a good 
statistical quality. Furthermore, the presented data stems from a well-resolved simulation (femax?? ~ 2), 
which is necessary as, e.g. , seco nd derivatives of the vorticity field will be considered. We also refer to the 
discussion in the literature^^ES regarding the importance of the spatial resolution. 



V. RESULTS 



Before coming to the investigation of the conditional vorticity budget, we show some benchmark results. 
Details on the CVE and Fourier filtering are summed up in tables |lT] and |III[ For example, it can be seen 
that with about 2.3% retained coefficients still about 99% of the enstrophy is contained in the coherent 
contribution. Due to the high cutoff wavenumber the Fourier filtering performs comparable. Furthermore, 
CVE represents the minimum and maximum values of the vorticity components to a better extent than the 
Fourier filtering. In figure [T] volume visualizations of the total, coherent, incoherent, low pass filtered and 
high pass filtered contributions of the vorticity are shown. It can be seen that both the coherent and low pass 
filtered contributions represent the global structure of the total vorticity field to a good extent, differences 
are only visible in the details. When comparing incoherent and high pass filtered contributions, differences 
become more pronounced. While the incoherent contributions appear very noisy and small in amplitude, 
the high pass filtered contributions have a much larger amplitude and appear more structured. The color 
scales for the visualization of the incoherent and high pass filtered fields have been adjusted to account for 
that issue. A zoom into the fields, however, also indicates differences in the low pass filtered contribution 
compared to the total field, as presented in figure [2J Although the shape of the vortex structures is captured 




Figure 2. Same visualization as in figure IT] but in a close-up view. Both the coherent and low pass filtered 
contributions represent the structure of the vorticity field rather well, the low pass filtered contribution tends to 
underestimate, for example, the core of coherent vortex structures. These contributions are contained in the high 
pass filtered contribution. 
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Figure 3. Left: PDFs of the magnitude of the total vorticity, the coherent part and the incoherent part of the wavelet- 
decomposed fields. Right: PDFs of the total, low pass filtered and high pass filtered Fourier-decomposed fields. The 
PDFs are highly non-Gaussian with stretched exponential tails in the case of the total and coherent vorticity. The 
incoherent part displays a much lower amplitude. Inset: Close-up of the PDF of the incoherent/high pass filtered 
vorticity. The nearly exponential/stretched exponential tail of the PDF indicates that also the incoherent/high pass 
filtered part is non-Gaussian. Note that the PDFs have been normalized by the standard deviation of the total 
vorticity a. The magnitude PDF corresponding to a Gaussian distributed vorticity field with standard deviation a 
is shown for reference. 
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Figure 4. Longitudinal vorticity autocorrelation functions {ujx{x) ujx{x + rex)) / {(-dx{x)^) for the total, coherent/low 
pass filtered and incoherent/high pass filtered vorticity (left: wavelet decomposition, right: Fourier decomposition). 
While the total and coherent/low pass correlations almost coincide, the autocorrelation function of the incoher- 
ent/high pass filtered part is rapidly decaying and oscillating. The incoherent part of the wavelet-decomposed 
fields is longer correlated than the high pass filtered part of the Fourier-decomposed fields. It can be seen that the 
incoherent/high pass filtered fields are much shorter correlated than the coherent/low pass filtered contributions. 
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Figure 5. Enstrophy spectra for the total, coherent and incoherent contribution as well as the large-scale and 
small-scale contributions. The coherent spectrum matches the total one over a broad range of scales, small-scale 
deviations, however, are visible. The Fourier filter (vertical line at fcc = 91) clearly separates large-scale and small- 
scale contributions. The high-frequency contributions of the wavelet spectra in the coherent and incoherent flow are 
due to a change of the basis functions, i.e. from trigonometric polynomials (Fourier basis) to Coiflet 30 (wavelet 
basis), which are not compactly supported in spectral space. 



quite well, it can, for example, be seen that the amplitude of the vortex cores sometimes is underestimated 
in the low pass filtered contribution. 

This observation can be made more quantitative by investigating the PDF of the magnitude of the vorticity 
which is presented in figure [3| This figure shows that the PDF of the coherent part of the vorticity yields a 
PDF almost indistinguishable from the total PDF. The PDF of the incoherent part has a largely reduced 
variance and displays a nearly exponential decay consistent with previous findingsISHiSI, Compared to that, 
the Fourier filter does not perform as well because the PDF of the low pass filtered part of the vorticity 
field is clearly distinguishable from the total PDPi^l, Furthermore, the high pass filtered component of 
the vorticity field has a much larger variance compared to the PDF of the incoherent part and displays a 
stretched exponential shape. 
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Figure 6. Conditional balance of the conditional averages related to vortex stretching and diffusion of vorticity. 
Vortex stretching is positively correlated with vorticity, whereas the diffusion of vorticity is negatively correlated. As 
expected, the sum of both terms cancels, indicating that the conditional balance ( 12 I holds. The gray lines indicate 



the functional form of the conditional averages expected when assuming statistical independence of the rate-of-strain 
tensor and the vorticity magnitude and a corresponding balance of the diffusive term. All conditional averages 
are normalized by {e^)/a consisting of the enstrophy dissipation and the standard deviation of the vorticity field, 
respectively. 
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Figure 7. Coherent and incoherent parts of the diffusive term compared to the total one (left: wavelet decomposition, 
right: Fourier decomposition). In the case of the wavelet decomposition the coherent part contributes the most, 
however, the incoherent part is small, but non- vanishing. In the case of the Fourier-decomposed fields the low pass 
filtered and high pass filtered contributions do not separate that well. As a benchmark, the sum of both contributions 
is shown to add up to the total diffusive term, which has been calculated for reference. 



To quantify the two-point statistics of the field, figure |4] shows the longitudinal correlation function of 
the different contributions. It can be seen that both the coherent and low pass filtered contributions of the 
vorticity field almost coincide with the correlation function of the total field. Differences, however, become 
apparent when comparing the incoherent and high pass filtered contributions. Although both correlation 
functions decay rapidly (compared to the correlation function of the total field), a long-ranging oscillation 
of the autocorrelation function of the incoherent field can be observed. These oscillations may also be seen 
in the visualizations (figures [1] and [2]) in form of a very fine-scaled structure of the incoherent field. 

These observations are also supported by studying the enstrophy spectra of the different contributions, 
c.f. figure [5] The enstrophy spectra of the total and coherent flows perfectly superimpose all along the 
inertial range. In the dissipative range, for wavenumbers larger than k > 60, we observe a departure, i.e.. 
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Figure 8. Coherent and incoherent contributions to the conditional average related to the vortex stretching term 
(left: wavelet decomposition, right: Fourier decomposition) . For the wavelet-decomposed fields the coherent-coherent 
contribution is dominant and almost identical to the total term. This indicates that vortex stretching is predominantly 
caused by both the coherent vorticity and the rate-of-strain tensor induced by the coherent vorticity. For the Fourier- 
decomposed fields the low pass-low pass contribution deviates significantly from the total one. 
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Figure 9. Contributions to the conditional vortex stretching term with at least one incoherent quantity (left: wavelet 
decomposition, right: Fourier decomposition). The conditionally averaged coherent part of the rate-of-strain tensor 
times the incoherent part of the vorticity contributes positively to the conditionally averaged vortex stretching term. 
In contrast to that, the term involving the incoherent rate-of-strain tensor times the coherent vorticity tends to 
deplete the vorticity. The term with both quantities incoherent seems to be negligible compared to the remaining 
terms. Note that the amplitude is significantly lower in the case of the wavelet-decomposed fields. 



the enstrophy spectrum of the coherent flows decays faster than the one of the total flow. In contrary, the 
enstrophy spectrum of the incoherent flow has a much weaker ampHtude and exhibits a slope close to fc"*. 
After reaching its maximum value at fc = 90, it rapidly decays. For reason of comparison, we also plotted the 
cutoff wavenumber kc = 91 corresponding to the black vertical line; this line divides the enstrophy spectrum 
of the total flow into large-scale contributions k < kc and small-scale contributions k > kc- According to the 
Wiener-Khinchin theorem, the spectra correspond to the Fourier transform of the autocorrelation functions. 
Hence, the oscillations observed in the longitudinal vorticity autocorrelation functions (see fig|4]) for both, 
the incoherent and the small-scale contribution are related to the maximum values in the corresponding 
spectra, i.e., the wavenumber of the oscillations is given by /cmax — (27r)/Zr where Ir denotes the wavelength 
of the oscillations and /cmax the maximum in the spectrum. 

We now come to the conditional vorticity budget and start with an investigation of the functional form of 
s{U) and d(0) presented in figure |6] It can be seen that the vortex stretching term is positively correlated 
with the vorticity, whereas the diffusive term is negatively correlated. This is physically quite intuitive as 
it mirrors the fact that the vortex stretching term tends to amplify vorticity, while the dissipative term 
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Z 


%Z 


min ijJi 


max uii % of retained coefF. 


total: 


51.6 


100 


-118 


123 100 


coherent: 


51.0 


98.8 


-118 


122 2.31 


incoherent : 


0.643 


1.25 


-10.7 


10.5 97.7 



Table 11. Enstrophy and percentage of retained coefficients of the total, coherent and incoherent contributions ob- 
tained using CVE. Mean value of the threshold e — 9.74. toe denotes the components of vorticity. All values are 
averaged over the fields of the ensemble. 





Z %Z 


min LJi 


max oji 


% of retained coeff. 


total: 


51.6 100 


-118 


123 


100 


low pass filtered: 


50.1 97.2 


-96.9 


98.1 


2.41 


high pass filtered: 


1.47 2.85 


-47.7 


49.5 


97.6 



Table III. Enstrophy and percentage of retained coefficients of the total, low pass and high pass filtered contributions 

obtained using Fourier filtering with cutoff wavenumber kc — 91, which was calculated as kc = {^N ■ Cr) ~ 91, 
where Cr = 2.31% is the percentage of retained coefficients for the CVE. All values are averaged over the fields of 
the ensemble. 
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0.978 


2.75 
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-5.85 
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0.913 


8.70- 10-2 


Fourier: 




0.931 


8.35 


10- 


_2 


-1.48 


10" 


-2 


4.62 


10" 


-4 




0.764 


0.236 



Table IV. Relative contributions of the different decomposed terms of the budget equation (121 to the average 
enstrophy budget. 



depletes vorticity. The fact that the sum of both averages nearly identically vanishes represents a posteriori 



justification for the approximation leading to the relation (12 1. In the same figure the functions expected 
for the case where the rate-of-strain tensor is assumed statistically independent of the vorticity and the 
corresponding diffusive term balances this term are shown for comparison. The slope of these linear functions 



is obtained such that these functions yield the correct ordinary enstrophy budget (13). The difference 
compared to the functions obtained from the DNS demonstrates, as expected, that pronounced correlations 
between the fields of the rate-of-strain tensor, the Laplacian of the vorticity, and the vorticity, respectively, 
exist. 

To now quantify the contributions of the coherent structures, we start with investigating the diffusive 
term, which is presented in figure [7J It is observed that this term is almost fully represented by the coherent 
part of the field, while the incoherent contribution appears significantly smaller. As the Laplacian of a field 
enhances its small-scale features, this demonstrates that the CVE captures these features especially well. 
This is also supplemented by the observation made for the Fourier decomposition. It can be seen that the 
high pass filtered component is smaller, but not negligible for low magnitudes of vorticity, which means that 
both fields contribute to the enstrophy dissipation. To make this more quantitative, these contributions have 



been calculated and are presented in table IV It can be seen there that the low pass filtered component 
contributes only about 76% to the total dissipation of enstrophy compared to 91% in the case of the coherent 
vorticity. 

Similar observations can be made for the terms related to the conditional vortex stretching term, which 
are shown in figure [Sj Also for this term the coherent contribution matches almost perfectly the total contri- 



bution, about 98% of the enstrophy production is contained within this term (see table IV ) . The remaining 
terms are strongly reduced in amplitude, still an investigation of their comparably small contributions is 
interesting, as can be seen in figure [9] It becomes apparent from this figure that the interaction of the 
coherent part of the rate-of-strain field with the incoherent part of the vorticity field is positively correlated 
with the vorticity, i.e., a positive contribution to the average enstrophy budget originates from this term. 
An interesting interpretation of this observation is that the rate-of-strain field produced by the coherent 
vortex structures is able to produce additional coherent vortex structures; in a sense coherent structures 
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breed coherent structures. In contrast to that, the interaction of the rate-of-strain field induced by the 
incoherent vorticity with the coherent vorticity has a depleting effect, such that it can be concluded that 
the incoherent rate-of-strain field destroys coherent vortex structures and hence has a dissipative character. 
The contribution of the incoherent rate-of-strain tensor times the incoherent vorticity field is negligible in 
view of its very low amplitude. 

For the Fourier filtering it can be seen that the low pass filtered contribution of the rate-of-strain field times 
the low pass filtered contribution of the vorticity field does not fully represent the total vortex stretching 



term. However, still about 93% of enstrophy production is contained within this term (see table IV) 



Although distinct differences compared to the CVE are apparent, table IV shows that these differences 
almost vanish in the average. This exemplifies that more detailed insights can be obtained by studying 
conditional averages instead of ordinary averages. For the cross-terms similar observations can be made as 
in the case of the wavelet analysis. The contribution of the low pass filtered rate-of-strain field and the high 
pass filtered vorticity field is positively correlated with the vorticity, whereas the opposite is observed for the 
case of the high pass filtered rate-of-strain tensor contribution and the low pass filtered vorticity field. Apart 
from the comparison to the wavelet filtered data, this is a physically interesting observation on its own as 
it quantifies the interaction of large- and small-scale contributions within the flow fields. The contributions 
involving high pass filtered fields also remain small in the case of the Fourier-filtered data. Still, they are 
significantly larger than the corresponding terms of the wavelet decomposition; the scale in figure [9] differs 
by more than an order of magnitude between the wavelet and Fourier decomposition. 



VI. CONCLUSION 

To summarize, we presented a detailed analysis of the conditional vorticity budget in terms of coherent 
vorticity. For this purpose we made use of the CVE technique to separate the noisy incoherent contributions 
of the vorticity field from the coherent ones. It was shown, in accordance with previous results, that CVE 
yields an excellent representation of the total flow using a strongly reduced number of degrees of freedom. 
This is particularly interesting as the conditional budget of vortex stretching and vorticity diffusion represents 
a dynamical rather than a purely kinematic relation. To further quantify the performance of the wavelet 
filtering method, we have performed a comparison to an ideal Fourier filter, where the fields have been low 
and high pass filtered. 

Although visualizations did not show strong discrepancies between Fourier and wavelet filtering, the in- 
vestigation of the conditional averages revealed more pronounced differences. It has been shown that most 
of the enstrophy production can be accounted to the coherent vorticity and the correspondingly induced 
rate-of-strain field. Interestingly, we have found that the incoherent rate-of-strain field tends to deplete 
vorticity, i.e., it tends to destroy coherent structures, while the rate-of-strain field induced by the coherent 
vorticity contributes positively. However, these contributions are small compared to the coherent-coherent 
contribution. In this sense the coherent structures are able to maintain or even amplify themselves, whereas 
the incoherent contributions tend to have a dissipative effect. Our analysis hence also quantifies the interac- 
tion of coherent and incoherent contributions to the flow. Consistent results have been found for the Fourier 
filtering. It has been shown, however, that this kind of filtering does not separate coherent contributions as 
well as wavelet filtering does, which is in agreement with previous studies. 

These findings motivate to further develop CVS, where neglecting the incoherent flow contributions is 
assumed to be sufficient to model turbulent dissipation. With respect to the performance of this method, 
it should be noted that Coherent Vorticity Simulations using adaptive wavelet bases require, however, the 
use of a safety zone to account for the generation of wavelet coefficients in scale and space due to the 
nonlinear flow dynamics. This increases the number of retained coefficients in actual simulations by a factor 
2-8 depending on the choice of the safety zone. For a discussion on possible choices and their infiuence on 
the flow statistics we refer the reader to a recent work^^. The Lagrangian particle-wavelet method^^ would 
allow to perform CVS without safety zone. Note that in simulations using Fourier spectral methods the 
number of modes has also to be increased by a certain factor in each spatial direction depending on the used 
dealiasing technique to remove high wavenumber modes for computing the nonlinear term. 
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